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A Sierpinski number ID page 420] and |[l], is an odd positive integer, k, such that no positive integer 
in this infinite list is prime: 

fe2' + l,A:22+l,A:2^ + l,...,A:2" + l,.... (1) 

A Riesel number 1 1 1 is similar to a Sierpinski number, with — 1 replacing +1 in the above infinite 
list. Such a number is an odd positive integer, k, so that no positive integer in this infinite list is prime: 

A:2'-l,fe22-l,fe2^-l,...,A:2"-l,.... 

A cover, for such a A;, is a finite list of positive integers such that each integer, j, in the appropriate 
infinite list, has a factor, d, in the cover, with \ <d < j. 

Given a k and its cover, ACL2 is used to systematically verify that each integer, in the appropriate 
infinite hst, has a smaller factor in the cover. 

1 Introduction 

Sierpinski and Riesel numbers are not easy to find. To disqualify an odd positive integer as a Sierpiiiski 
number or a Riesel number, one need only locate a prime in the appropriate infinite list. With four 
exceptions, k = 47, 103, 143, 197, all of the first 100 odd positive integers, 1 < ^ < 199, are disqualified 
as Sierpinski numbers by finding at least one prime in the first eight elements of the infinite list |[3l : 

kl^ + \,kl^ + \,kl^ + 1, . . . ,k2^ + \. 

Both k = 103 and k = 197 are eliminated by finding a prime in the list no later than k2^^ + 1 |3], leaving 

47 and 143 as the only possible Sierpinski numbers less than 200. It turns out that 143 • 2^^ + 1 and 
47.2583 

are prime 131, eliminating them. Thus, there are no Sierpiiiski numbers in the range I <k< 
199. The situation is similar for Riesel numbers. 
In 1960, W. Sierpinski |7| proved, for 

;t = 15511380746462593381, 

every member in the infinite list, given by ([T|l, is divisible by one of the prime factors of the first six 
Fermat numbers. For nonnegative integer, n, the Fermat number, F„, is given by 

F„=2^ +1. 

The first five Fermat numbers are prime and is the product of two primes: 

Fo = 3,Fi = 5,^2 = 17,^3 = 257,^4 = 65537, 
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F5 = 4294967297 = 641 • 6700417. 

Thus (3 5 17 257 641 65537 6700417) is a cover for k = 15511380746462593381, showing it to be a 
Sierpiriski number. Sierpiriski's original proof is described in [8, page 374] and [2J. 

In 1962, J. Selfridge (unpublished) proved that 78557 is a Sierpiriski number by showing that 

(3 5 7 13 19 37 73) 

is a cover f6l. Later, in 1967, Selfridge and Sierpiriski conjectured that 78557 is the smallest Sierpin- 
ski number [6|. The distributed computing project Seventeen or Bust [6J is devoted to proving this 
conjecture, disqualifying every k < 78557, by finding an n that makes k-2" + 1 prime. For example, 
19249 •2i3"i^5^^ + 1, a 3918990-digit prime, eliminated 19249 Q. When this project started in 2002, all 
but 17 values of k had already been disqualified. Currently six values of k remain to be eliminated. 

Earlier, in 1956, but less well known than Sierpiriski's work, H. Riesel [5 1 showed 509203 is a Riesel 
number with cover (3 5 7 13 17 241). It is possible for the same odd positive integer to be both a 
Sierpiriski number and a Riesel number. An example LU is A: = 143665583045350793098657. 

2 Covers Into ACL2 Proofs 

Given an odd positive integer, k, with a Sierpiriski cover, here is the process used to verify that ^ is a 
Sierpiriski number. There is a similar process for verfying Riesel numbers from their covers. 

1 . For each d in ^, find positive integer bd and nonnegative integer so that for every nonnegative 
integer /, J is a factor of k ■ 2^''''^'^'' + 1. 

In practice, every J in is an odd prime smaller than k. 

(a) Search for positive integer b such that J is a factor of 2'' — 1. Since d is an odd prime, it turns 
out that such a b will always exisfl among 1,2,. — 1. Let b^ be the firsjl such b. 

(b) Search for nonnegative integer c such that <i is a factor of /: • 2' + 1. If such a c exists, then 
one exists among 0, 1 , . . . , Zj^ — 1 . Let Cd be the firsH such c, if it exists. 

(c) Assuming Cd exists, use induction on /, to prove that for every nonnegative integer /, (i is a 
factor ofk-l^'> '+'" + 1. 

The base case, when / = 0, follows from [Tb] above. 

The induction step, going from / = j to / = 7 ' + 1, follows from [Tal above: 

;t2^<i(y+i)+crf ^ y ^ [kl^^dj+^d . (2''^ _ 1)] + [k2''dJ+'-'' + 1] (2) 

By [Tal d is a factor of the left summand of dD and J is a factor of the right summand by the 
induction hypothesis. 

2. For each positive integer n, find d in'^ and nonnegative integer / so that n = bd -1 + Cd- If such d 
and / exist, then, by [TB is a factor of k ■ l^d-'+'^d + 1 = ;t • 2" + 1. 

To ensure that such d and / exist for every positive n, only a finite number of cases need be 
considered: Let I be the least common multiple of all the bd's found for the cf's in '10. Check for 
each 

«G{0,1,2,...,£-1}, 

'For the mathematically literate: The well-known Fermat's Little Theorem ensures the claimed existence. 
^Thus, being mathematically precise, i)^ is just the order of 2 in the multiplicative group of the integers modulo d. 
^If d does not divide k, then 2'^'' = —[\/k) (modd), so is the discrete logarithm, base 2, of — (1/A:) in the integers 
modulo d. 
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that there always is a J in that satisfies the equation 

mod in,b4) =q. 

This process has not been formally verified in ACL2. For example, we don't bother to check that every 
member of is an odd prime. Instead, for each individual k and ACL2 events are generated that 
would prove is a Sierpinski number, if all the events succeed. If some of the events fail, then, as 
usual when using ACL2, further study of the failure is required, in the hope of taking corrective action. 
The generation of these events is controlled by the macros verif y-sierpinski and verif y-riesel. 
These macros take three arguments: the name of a witness function that will find a factor for a given 
k2" ± 1, the number k that is a Sierpinski or Riesel number, and the cover ^ for k. The macros then 
generate the proof, following the plan outlined in this section. 

For each d in'j^, and q from la and lb, are computed. They are needed to define the witness 
function and to state the theorems mentioned in Ic, which are then proved. For example, the proof that 
78557 is a Sierpinski number defines this witness function: 

(DEFUN WITNESS (N) 
(IF (INTEGERP N) 
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0)) 

The rightmost numbers, in this definition, form the cover, the corresponding b/s are the leftmost num- 
bers, and the middle numbers are the q's. So ^ = (3 5 7 13 19 37 73), bjs = 9, and C73 = 3. 
The theorem, from Ic, for J = 73 is 

(DEFTHM WITNESS-LEMMA-73 
(IMPLIES (AND (INTEGERP N) 
(>= N 0)) 
(DIVIDES 73 

(+ 1 

(* 78557 
(EXPT 2 

(+ 3 

(* 9 N))))))) 

: HINTS . . . ) 

Four properties are proved about the witness function, establishing 78557 is a Sierpinski number: 

(DEFTHM WITNESS-NATP 
(AND (INTEGERP (WITNESS N)) 

(<= (WITNESS N))) 
: HINTS . . . ) 
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(DEFTHM WITNESS-GT-1 
(IMPLIES (INTEGERP N) 

(< 1 (WITNESS N))) 
: HINTS . . . ) 

(DEFTHM WITNESS-LT-SIERPINSKI 
(IMPLIES (AND (INTEGERP N) 
«= N)) 
(< (WITNESS N) 

(+ 1 (* 78557 (EXPT 2 N)))))) 

(DEFTHM WITNESS-DIVIDES-SIERPINSKI-SEQUENCE 
(IMPLIES (AND (INTEGERP N) 
(<= N)) 
(DIVIDES (WITNESS N) 

(+ 1 (* 78557 (EXPT 2 N))))) 

:HINTS . . .) 

As suggested above in 2, these properties can be proved by showing every integer is "covered" by one of 
the cases given in the COND-expression used in the definition of the witness function. 

(DEFTHM WITNESS-COVER-ALL-CASES 
(IMPLIES (INTEGERP N) 

(OR (EQUAL (MOD N 2) 0) 

(EQUAL (MOD N 4) 1) 

(EQUAL (MOD N 3) 1) 

(EQUAL (MOD N 12) 11) 

(EQUAL (MOD N 18) 15) 

(EQUAL (MOD N 36) 27) 

(EQUAL (MOD N 9) 3))) 
: RULE-CLASSES NIL 
: HINTS . . . ) 

To prove this, we first demonstrate that these cases are exhaustive when n is replaced by mod(?i,36) 
(where 36 is the least common multiple of all the moduli above). This can be checked, essentially, by 
computation. 

(DEFTHM WITNESS-COVER-ALL-CASES-MOD-36 
(IMPLIES (INTEGERP N) 

(OR (EQUAL (MOD (MOD N 36) 2) 0) 

(EQUAL (MOD (MOD N 36) 4) 1) 

(EQUAL (MOD (MOD N 36) 3) 1) 

(EQUAL (MOD (MOD N 36) 12) 11) 

(EQUAL (MOD (MOD N 36) 18) 15) 

(EQUAL (MOD (MOD N 36) 36) 27) 

(EQUAL (MOD (MOD N 36) 9) 3))) 
: RULE-CLASSES NIL 
: HINTS . . . ) 
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The actual modular equivalences that need to be proved depend on both the number 78557 and its cover. 
Although the theorem that is being proved is obviously true, there does not appear to be a way to prove it 
once and for all in ACL2, not even using encapsulate. Instead, a pair of theorems very much like the 
ones we have described needs to be proved from scratch for each different Sierpinski or Riesel number. 
As experienced ACL2 users, we are concerned that ACL2 will simply fail to prove this theorem for some 
combination of numbers and their covers. However, we have used these macros to generate the proof for 
each of the Sierpinski and Riesel numbers with covers listed in the appendix, and all of the proofs have 
gone through automatically. Note that the appendix essentiall>0 contains all the Sierpinski and Riesel 
numbers known to us. 

3 Numbers Without Covers 

There are odd positive integers, shown to be Sierpinski (or Riesel) numbers, that have no known covers. 
ACL2 proofs have been constructed for these numbers. 

For example id, A: = 4008735125781478102999926000625 is a Sierpinski number, but no (com- 
plete) cover is known. For all positive integer, n, if mod {n,4) ^ 2, then ^ • 2" + 1 has a factor among the 
members of (3 17 97 241 257 673). To show kisa Sierpinski number, a factor of ^ • 2" + 1 must be found 
for all positive integer, n, such that mod {n,4) = 2. Such a factor is constructed using these facts: 

• k = 44745755^^ is a fourth power 

• 4x'^ + \ = {2x^ + 2x+\)-{2x^-2x+\) 
Let / = 44745755, sok = Z^. Then 

k-2'^"+^ + i = 2^{i-rf + i 

= 4(/-2")Vl 

= [2(/-2")2 + 2(/-2") + l]-[2(/-2")2-2(/-2") + l] (3) 
The left factor of ^ algebraically reduces to show 

4004365 181040050 -22 L"/4J + 89491510 • 2L"/'^J + 1 

is a factor of • 2" + 1, whenever mod {n,4) = 2. 

A Riesel number, k, with no known cover, is given in Appendix A. In this example, k = a^ is a square 

and 

k-2^"-l = a^-2^"-l 
= {arf-l 
= {aT + I) ■ {a2" - I) 

shows how to factor ^ • 2*" — 1 when m is even and positive. A (partial) cover, listed in Appendix A, gives 
a constant factor for each k-2'" — I, when m is odd and positive. 

^Given a Sierpinski or Riesel number k and its cover infinitely many other examples can be constructed: Let P be the 
product of the numbers in "if and let i be a positive integer. Then k + 2-i-Pis also a Sierpinski or Riesel number with the same 
cover. 
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4 Conclusions 

Given a Sierpiriski or Riesel number, k, and its cover, we have described ACL2 macros that generate 
events verifying that each integer, in the appropriate infinite list, has a smaller factor in the cover. 

For the few known Sierpiriski or Riesel numbers with no known covers, hand-crafted ACL2 proofs 
have been constructed verifying that each integer, in the appropriate infinite list, has a smaller factor. 
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A Sierpinski and Riesel Numbers 

Numbers, k, verified with ACL2. 

Each k with a cover ^ is either mentioned in the References or claimed at various websites. Numbers 
k without known covers are from HI. 

Smallest known Sierpinski number 

k = 78557 = 17 • 4621, a product of two primes 

'r= (3 5 7 13 19 37 73) 

Smallest known prime Sierpinski number 

A: = 271 129 

'(^= (3 5 7 13 17 241) 
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More Sierpinski numbers 
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Smallest Sierpinski number found by Sierpinski 

yt = 15511380746462593381 

■^ = (3 5 17 257 641 65537 6700417) 

Smallest known Riesel number 

k = 509203 

■^ = (357 13 17241) 

More Riesel numbers 

k ^ 

762701 (3 5 7 13 17 241) 

777149 (3 5 7 13 19 37 73) 

790841 (3 5 7 13 19 37 73) 

992077 (3 5 7 13 17 241) 

Numbers both Sierpinski and Riesel 

'j^R indicates the Riesel number cover and indicates the Sierpinski number cover. 



k = 143665583045350793098657 

= (3 5 13 17 97 241 257) 
<^s = (3 7 11 19 31 37 61 73 109 151 331 1321) 



k = 47867742232066880047611079 
■Tr = (3 7 11 19 31 37 41 61 73 109 151 331) 
= (3 5 13 17 97 241 257) 
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k = 878503122374924101526292469 
-Tr = (3 7 13 19 37 73 97 109 241 257) 
= (3 5 11 17 31 41 61 151 331 61681) 



k = 3872639446526560168555701047 
-Tr = (3 7 13 19 37 73 97 109 241 673) 
-Ts = (3 5 11 17 31 41 61 151 331 61681) 



k = 623506356601958507977841221247 

= (3 7 13 19 37 73 97 109 241 673) 
"^s = (3 5 17 257 641 65537 6700417) 

Sierpinski numbers without cover 

k = 4008735125781478102999926000625 = 44745755^ 
(3 17 97 241 257 673) is partial cover for mod {n,4) / 2. 
4004365181040050 • 22L"/4J + 89491510 • 2L"/4J + i 
is a factor of A: • 2" + 1 , whenever mod (n, 4) = 2. 



A: = 734110615000775^ 

(3 17 257 641 65537 6700417) is partial cover for mod (riA) 7^ 2. 
10778367901 13632192906501201250 • 22L«/4J + 1468221230001550 • 2L«/4J + 1 
is a factor of A • 2" + 1, whenever mod (n,4) = 2. 

Riesel number without cover 

Let a = 3896845303873881175159314620808887046066972469809 and let k = a^. 

The hst 

(7 17 31 41 71 97 113 127 151 241 257 281 337 641 673 1321 14449 29191 65537 6700417) 

is partial cover for odd positve n. 

a ■ 2"/^ + 1 is a factor of A • 2" — 1, whenever n is positive and even. 



